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Abstract. We study an inverse problem for nonlinear elliptic 
equations modelled after the p-Laplacian. It is proved that the 
boundary values of a conductivity coefficient are uniquely deter- 
mined from boundary measurements given by a nonlinear Dirichlet- 
to-Neumann map. The result is constructive and local, and gives 
a method for determining the coefficient at a boundary point from 
measurements in a small neighborhood. The proofs work with the 
nonlinear equation directly instead of being based on linearization. 
In the complex valued case we employ complex geometrical optics 
type solutions based on p-harmonic exponentials, while for the real 
case we use p-harmonic functions first introduced by Wolff. 



1. Introduction 

This article concerns inverse boundary value problems for nonlinear 
elliptic equations. In the case where the underlying equation is linear, a 
standard example is the inverse problem of Calderon [13]. This problem 
is related to Electrical Impedance Tomography, a method proposed for 
medical and industrial imaging, where the objective is to determine 
the electrical conductivity of a medium by making voltage to current 
measurements on its boundary. There is an extensive theory concerning 
inverse boundary value problems for linear elliptic equations. We refer 
to [H] for a recent survey. 

We recall the mathematical statement of the Calderon problem. Let 
Q C R n , n > 2, be a bounded open set with C 1 boundary. We consider 
fl as a medium which conducts electricity, with conductivity given by 
a positive function 7 e L°°(Q). Assuming that there are no sources or 
sinks of current in O, a voltage / on the boundary induces a potential 
u in the domain which (by Ohm's law) solves the Dirichlet problem 

div(7(x)VM) =0 in ft, 
u = f on dfl. 
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The boundary measurement corresponding to /, denoted by A 7 /, is 
the current at the boundary given by 



We assume that one can prescribe many different voltages / on the 
boundary, and then measure the corresponding boundary currents A 7 /. 
The map A 7 is called the Dirichlet-to-Neumann map (DN map for 
short). Using a suitable weak formulation, A 7 becomes a bounded 
linear operator from W 1 ^ 2,2 (d^l) to W~ l / 2,2 (dQ) (where W s ' p denotes 
the LP based Sobolev space with smoothness index s). The Calderon 
problem asks to determine the conductivity function 7 from knowledge 
of the operator A 7 . 

There are several aspects of the Calderon problem that have been 
studied. We mention the following particular questions: 

1. Boundary uniqueness. If A 7l = A 72 , then 7i|an = 72)9^- 

2. Interior uniqueness. If A 7l = A 72 , then 71 = 72 in Q. 

3. Stability. If A 7l and A 72 are close, then 71 and 72 are close. 

4. Reconstruction. Algorithm for determining 7 from A 7 . 

5. Partial data. If A 7l |r = A 72 |r for some T C dQ, then 71 = 72. 

As a rule, rather precise results for the above questions are available 
in the case n = 2 [TTJ, [8], [9], [TJ, [17]. Also in dimensions n > 3 
many results have been obtained [TTJ, [39], [TJ, [31], [26], however sharp 
conditions such as optimal regularity of the conductivity are in general 
not known. We refer again to the survey (TTJ for more details. 

In contrast with the linear case, less is known about variants of the 
Calderon problem for nonlinear elliptic equations. In this paper we 
consider a particular nonlinear model based on the p-Laplace operator 



The corresponding p-Laplace equation A p u = 0, whose solutions are 
called p-harmonic functions, is a prototypical nonlinear equation in di- 
vergence form. It arises as the Euler-Lagrange equation for minimizing 
the p-Dirichlet energy 



over all u G W l,p (Q) with fixed boundary values. For more details on 
p-Laplace type operators we refer to the book [15] and lecture notes 
[Lf] . [30J. Applications in fluid mechanics, plastic moulding, and image 
processing are discussed in [TJ, [5], [28] . 



A p u = div(|Vu| p - 2 



Vw), 1 < p < 00. 
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Given a bounded open set Q C M. n , n > 2, having C 1 boundary, and 
for a positive function 7 e L°°(f2), we consider the Dirichlet problem 

f div( 7 (x)|Vu| p " 2 Vu) = infi, 
1 it = f on <9f2. 

This is also called the weighted p-Laplace equation [15] . The nonlinear 
DN map is formally defined by 

A 7 : / h-> 7 (x)|Vm| p ~ 2 Vm • i/|an 

where it is the unique solution with boundary values / and v is the 
outer unit normal to dQ. The precise definition of the DN map is 
given in Appendix lAl 

Since the equation is nonlinear, one needs to make a distinction 
between real and complex valued solutions. We denote by and A^ 
the DN maps acting on real and complex boundary values, respectively. 
Our main theorem is the following boundary uniqueness result. 

Theorem 1.1. Let Q C W 1 , n > 2, be a bounded open set having 
C l boundary, and let 71,72 be positive continuous functions on Q. If 
A^ = A*, ^en7i|en = 72|an. 

The proof is constructive and local in the sense that we construct 
a sequence of explicit functions fu on the boundary, supported in an 
arbitrarily small neighborhood of a boundary point Xo, such that 

lim / k 1 (f M )f M dS = 7(x ). 

In fact, this result is true if 7 is in L°°(0) and continuous near x - The 
map A^ determines A^ trivially, and we obtain the following conse- 
quence for the complex valued case. 

Theorem 1.2. Let Q C M. n , n > 2, be a bounded open set having 
C 1 boundary, and let 71,72 be positive continuous functions on Q. If 
A 7i = A 72? then 7i I an = l2\an- 

The reason for stating the theorems separately is that we actually 
first prove Theorem 11.21 by employing solutions based on p-harmonic 
complex exponentials, and then establish Theorem II .11 by isolating the 
properties of complex exponentials needed for the proof and by making 
use of suitable real valued p-harmonic functions. 

There are many works concerning boundary determination in the 
linear case where p = 2. It was proved in [27] that the Taylor series 
at a boundary point of a smooth conductivity in a smooth domain is 
determined by the DN map. Another proof, based on pseudodiffer- 
ential calculus and valid in many situations, was given in [ID]. For 
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nonsmooth domains and conductivities there are boundary uniqueness 
results based on singular solutions [2], [3] and direct methods involving 
explicit oscillating boundary values [IT] . [25] . [32] . [33] . 

Let us describe earlier results on nonlinear variants of the Calderon 
problem. One can consider the Dirichlet problem involving a more 
general nonlinear conductivity a = a(x, z, q), 

div(a(x, u, Vw)Vw) =0 in Q, 
u = f on dfl. 

Examples include 

1. (linear case) a(x, z, q) = j(x) for a positive 7 G C 2 (f2), 

2. (nonlinearity depending on x and u) a(x, z, q) = 7(2, z) where 
7 G C 2 ((] x 1) is positive, 

3. (derivative nonlinearity) a(x, z, q) = j(x, q) for suitable 7. 

In all these cases, the equation is uniformly elliptic and there is a unique 
solution in a suitable Sobolev space for any suitable boundary value /. 
We can then define the nonlinear DN map formally by 



A 7 : / \-> a(x, u, Vw) Vw • v\dti- 



As explained above, many results are known for the linear case. Also 
the case where the nonlinearity depends on x and u is well understood. 
In fact, it was shown by Sun [M] that this case reduces to the linear 
theory (the linearization idea is due to [H]). If 1 < p < 00, if z G E, 
and if 7 z (x) = 7(2;, z), then one can prove that for any / G C 2,a (dQ) 

lim A^-H/)-A 7 (z) 

i^O t 7 

in W l - l/p ' p {dtt). As a consequence, if A 7l = A 72 for two conductivities 
7j = 7j(sc, m), then A 7 ^ = A 7 | for all z G M. and the interior uniqueness 
result in the linear case shows that 71 = 72 everywhere. Related results 
for other equations with nonlinearity depending on x and u appear in 

HB|, HSI, |2U, 122], 

For derivative nonlinearities it is possible to obtain some information 
from linearizations, see [IB] . [23] . [3B] for conductivity type equations 
and [20] , [2H] , [32] for related equations. These results are still based on 
linearizing the nonlinear DN map and applying the known uniqueness 
results for the linear case, and they apply to derivative nonlinearities 
of certain form. To deal with stronger nonlinearities one could hope 
for a method which works with the nonlinear equation directly. In this 
paper we introduce such a method, at least for the purposes of proving 
boundary determination results for p-Laplace type equations. 
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In the linear case most uniqueness results are based on special solu- 
tions of the conductivity equation, so called complex geometrical optics 
solutions, which look like harmonic exponentials e px where p G C n 
and p ■ p = 0. The first important observation is that also for the 
p-Laplacian there exist special complex solutions of the form e p ' x (see 
Lemma 12. ip . We prove Theorem 11.21 by showing that one can per- 
turb these p-harmonic complex exponentials to become solutions of 
the equation involving 7, concentrating near a boundary point. The 
proof is similar to the arguments of [UJ, [T2] in the linear case and is 
actually not that difficult, making use of basic facts such as wellposed- 
ness for the Dirichlet problem, inequalities for pth powers of vectors, 
and Hardy's inequality. We then show Theorem 11.11 by replacing the 
p-harmonic complex exponentials with certain real valued p-harmonic 
functions, introduced by Wolff [42J, having similar properties as expo- 
nentials which allow the proof to go through. 

It is an interesting question whether one can make progress on other 
aspects of inverse problems for p-Laplace type equations (such as inte- 
rior uniqueness, stability, reconstruction, partial data) besides bound- 
ary uniqueness. Our results seem to suggest the possibility of a com- 
plex geometrical optics construction based on p-harmonic complex ex- 
ponentials, or a corresponding construction in the real case using the 
functions of Wolff. One also expects the cases n = 2 and n > 3 to 
be different. (In particular, when n = 2 the p-harmonic equation is 
related to quasiregular mappings [6] and also unique continuation for 
p-Laplace is known when n = 2 [10] but it is not known for n > 3.) 
Another interesting direction would be to study more general equations 
modelled after the p-Laplacian, or to see if methods of this type are 
available for other physically relevant nonlinear equations. This paper 
is mainly intended to highlight a particular strongly nonlinear model 
for which Calderon type problems can be studied, and to give a first 
result in this setting. 

This paper is organized as follows. Section [TJ is the introduction. 
The complex valued case is considered in Section [2] where Theorem 11.21 
is proved, while Section [3] discusses the real valued case and proves 
Theorem 11.11 For the sake of completeness, there is an appendix con- 
taining standard material on inequalities for pth powers of vectors and 
on wellposedness and the Dirichlet-to-Neumann map for the equations 
considered in this paper. In the appendix we also make the observa- 
tion that linearization of the DN map at constants does not give useful 
information of the conductivity, thus partly justifying the nonlinear 
methods used to prove the theorems. 
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Notation. If z, w G C n we write z-w = Y^j=i z j w j f° r the dot product 
and \z\ = (z-z) 1 / 2 = (|Re(2)| 2 + |Im(2;)| 2 ) 1 / 2 for the norm. If 1 < p < oo 
the conjugate exponent is denoted by p', so that 1/p + 1/p' = 1. We 
write W s,p for the standard Sobolev spaces and Wq' p for the closure of 
smooth compactly supported functions in W l,p . The notation A < B 
means that A < CB for some constant C > which is independent 
of asymptotic parameters (it typically only depends on n, p, and some 
choices of test functions). Similarly, A ~ B means that ±A < B < CA 
for some constant C > 0. 

Acknowledgements. M.S. is supported in part by the Academy of 
Finland, and is grateful to the Department of Mathematics and Statis- 
tics of the University of Jyvaskyla where part of this work was carried 
out. X.Z. is supported by the Academy of Finland. 

2. Complex case 

We now prove Theorem 11.21 concerning boundary determination from 
the DN map with complex boundary values. To convey the main idea 
without unnecessary complications, we will consider the case where Q 
is a bounded open set in M. n with C 1 boundary and xq is a point in 
dfl such that dfl is flat near xq. (The non-flat case is covered when 
proving Theorem 11.11 in Section |3j) By a translation and rotation, we 
may assume that x = and fl PI -8(0, r) = {x G -8(0, r) ; x n > 0} for 
some small r > 0. 

If 7 = 7(2;) G C(Q) is a positive function, consider the Dirichlet 
problem 

f dw{-f{x)\Vu\ p - 2 Vu) =0 in Q, 
1 u = f on <9f2. 

The nonlinear DN map, acting on complex valued functions, is defined 
in the weak sense (see Appendix [A} by 

/ A^(f)gdS= [ -?\Vu\ p - 2 Vu-Vvdx, f^eW 1 *^) 
Jan Jn 

where v is any function in W 1,p (fl) with v — g G Wq ,p (Q). 

The main idea in the proof is to use a special solution to the non- 
linear equation with coefficient 7 frozen at the boundary point 0. This 
equation is just the p-Laplace equation, and the special solution is the 
following p-harmonic complex exponential. 

Lemma 2.1. Let h(x) = e p ' x where p = a + i/3 with a, (3 G MJ 1 . Then 
h satisfies divfl Vh\ p ~ 2 Vh) = iff (p - l)\a\ 2 = \/3\ 2 and a- (3 = 0. 



AN INVERSE PROBLEM FOR THE p-LAPLACIAN 



7 



Proof. Since Vh = pe p ' x , we have 

div(|V/i| p - 2 V/i) = dw{\p\ p - 2 e ip ~ 2)ax pe p - x ) 

= dw{\p\ p - 2 pe ip - 1)a - x+il3 - x ) 

= \P\ P ~ 2 P ■ {(P - l)a + if3) e {p ' 1)a - x+ip - x . 

Here p • ((p — l)a + = (p — 1)|«| 2 — |/3| 2 + «pa • which proves the 
result. □ 

We wish to convert the p-harmonic function e p ' x into an exact solu- 
tion of div(7| Vu| p ~ 2 Vu) = in Q which concentrates near the bound- 
ary point 0. To this end, define the function 

(2.1) u (x) = r) M {x)h N (x) 

where T]m(x) = r)(Mx), h^ix) = h(Nx) where M and N are large 
positive numbers, rj G C%°(M. n ) is a nonnegative cutoff function with 
7] = 1 for \x\ < 1/2 and rj — for \x\ > 1, and 

h(x) = e (i(3 - en) - x 

with P G M. n satisfying \(3\ 2 = p — 1 and /3 ■ e n = 0. We will choose 
N = N(M) so that M/N as Af oo. The idea is that with these 
choices, since solves the equation with 7 frozen at and since Uq is 
supported in the ball B(0, 1/M), wo becomes an approximate solution 
to the nonlinear equation in Q when M is large. Lemma 12.41 below 
gives a precise meaning to this statement. 

We obtain an exact solution u by solving the Dirichlet problem with 
boundary values u , 

( div( 7 (a;)|Vw| p - 2 W) = in Q, 
1 u = u on dVL. 

Let / = Uo\an- Then we have 

/ A 7 (f)fdS= [ j\Vu\ p - 2 Vu-Vu dx. 
Jon Jn 

We write this as 

(2.2) / A,(f)fdS= [ j\Vu \ p dx 
Jan Jn 

+ / 7(|Vu| p - 2 Vu - |Vn |^ 2 Vn ) ■ Vu dx. 
Jn 

Note that since / is an explicit function, the left hand side is determined 
by the nonlinear DN map. We will recover the value of 7 at by taking 
the limit of this identity as M — » 00. To analyze the limit, we need a 
simple lemma. 
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Lemma 2.2. Let ( e C™(B(Q, 1)) and let a > 0. Then as M -> oo 
M n_1 iV / C(Mx)e- pNx " dx -» - / C(»', 0) 
i^C(Mx)e- pJVl " cfe = 0(M 1_ "iV~ 1_a ). 



Proof. Follows from a direct computation. □ 
We compute the limit of the first term on the right hand side of 

Lemma 2.3. We have as M — )■ oo 

M n-i N i- P f 7 | V u | p rfa; ->• c p7 (0) 
in 

where c p = p~ L„_i ^(^', 0) p <ix / . VFe a/so /tai>e 

|Vw | p rfx = 0(M 1 - n A^ p - 1 ), 



|77 M V/ijvr^ = 0(M 1 - n A^ p - 1 ), 

/ |Vw - %/V^| p rfx = 0(M 1 ~ n N p ~ 1 (M/N) p ). 
Jn 

Proof. Since u = 77M^jV, we compute 

Vw = MVr?(M ■ )h N + r} M Vh N . 
Since V/iat = N(i(3 — e n )e N! - ll3 ~ en ^ x = N(i/3 — e n )h N , we have by Lemma 



\\MV V (M-)h N \\ p LP{n) = 0{M 1 ~ n N' 1 M p ), 
\\VM^h N \\ p LP{Q) = 0(M l - n N- l N p ). 

This shows the last three estimates since M/N = o(l) as M — > oo. 
For the first statement, we use the inequality (IA.6I) to conclude that 

/ {\Vu \ p -\N{i(5 - e n )n M h N \ p ) dx 
Jn 

<p [ ^V^M -)h N \{\Vu \ p ~ l + \N{tf3 - e n ) m h N \ p ' l )dx 



<p||MVr/(M-)^||L^)(||VM || p ; ( 1 f , ) + \\N{if3 - e n )r]Mh N \\ p -l n) ) 

rl-n j\rp-l 



0(M L - n N p - L (M/N)). 
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Using that \i(3 — e n \ 2 = p, we have by Lemma I2T21 

lim M n ^ 1 N 1 ^ p [ \Vu \ p dx= lim M n ~ 1 Np p ^ 2 [ \r] M h N \ p dx 

r](x',0) p dx'. 



P-2 
P 2 



The result follows by writing 7 = 7(0) + (7 — 7(0)) and by using the 
continuity of 7. □ 

We now move to the analysis of the second term on the right hand 
side of (12. 2p . Writing u = uq+ui, the next result shows that || Vwi||,LP(n) 
is asymptotically smaller than || VwoIIlp^)- This may be interpreted 
so that Mi is a small correction term which corrects the approximate 
solution Mo into an exact solution u. The important facts for the proof 
are that A p h = and that Uq is supported near the boundary which 
makes it possible to use Hardy's inequality: if 5(x) = dist(x, dQ) then 

\\v/S\\ LP{n) < C\\Vv\\ LP{n) , v e W^ P (Q). 
Lemma 2.4. As M — > 00 

/ \V Ul \ p dx = o(M 1 ~"iV p - 1 ). 
Jn 

Proof. We will prove that 

I=[ (|Vw| + \Vu \) p ~ 2 \V Ul \ 2 dx 
Jn 

(2.3) < o(M l ~ n N p ~ l ) + o(l) [\V Ul \ p dx. 

Jn 

To prove fl 2 . 3 j) . we start with 

i< [ i(\vu\ + \vu \y- 2 \v Ul \ 2 dx, 

Jn 



since 7 is positive on Q. Then we invoke the inequality (IA.8j) . Since 
u\ = u — uq € jy o 1,p (f2) and since u is a solution, we obtain that 



= - Re 



/ 7(|Vu| p " 2 Vu - |Vu | p " 2 Vu ) ■ (Vw - Vw ) dx 
Jn 

/ 7|Vu | p ~ 2 Vm ■ Vuidx 
Jn 
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The function w is supported in the ball B(0, 1/M). Consequently, 
writing 7 = 7(0) + (7 — 7(0)), we have 



/ < 



Vu \ p 2 V«n • Vwi dx 



u ■ \U\ 



+ I 17-7 
's(o,i/M)nn 

-h + h 



Integral J 2 is bounded by H7 - tC ) IU~(s(o,i/AOnn) II Vwollz/llVuiHiP, 
which implies by Lemma [2731 the continuity of 7 and Young's inequality 
that 

I 2 < o(M 1 - n A^ p - 1 ) + o(l) f \V Ul \ p dx. 

Jn 

Then we estimate integral I\ as follows. At this point it is convenient 
to replace V«o with ^AfV/ijv by writing 

Vm | p ~ 2 V'Uo • V«i dx = \rjM^h N \ p ~ 2 rjM^h N ■ Vwi dx 
n Jn 

+ / (|Vw | p ~ 2 Vw - \VM^h N \ p - 2 r] M Vh N ) ■ S/uidx. 
Jn 

Integrating by parts, we obtain that 



div(|?7 M V/iAr| p 2 r] M Vh N )uidx 



+ 



/ (|Vw r 2 Vw - \VM^h N \ p - 2 r] M Vh N ) ■ Vtii dx 
Jn 



In the first term on the right, we multiply and divide by 5 (the distance 
to the boundary) and use the Holder and Hardy inequalities so that 



div(|?7MV/iAr| p 2 r]M^h N )uidx 



< ||5div(|?7MV/lAr| P ?7MV/iAr)|| LP '||Vui|| L p. 

The second term on the right can be estimated by (1A.7I) . and we have 

/ (|Vu | p ~ 2 Vm - \VM^h N \ p - 2 r] M Vh N ) ■ Vwi dx 
Jn 

< / (|V«o| + \r] M Vh N \) p - 2 \Vu - ?7AfV/i7v||V«i| dx, 
Jn 
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which, by the Holder inequality, is bounded by 

(||Vm ||lp + ||^MV/lAr||z / p) P " 2 ||VM - ??M V/ijV || L p || V«i || LP 

= 0{{M l - n N p - l )*irM/N)\\Vux\\LP, 

when p > 2, and by 

||Vw - Vm^Hn I l^p 1 1| V«i || L p 

= 0((M 1 ^N v ~ 1 ) 1L r(M/N) p - 1 )\\\7u 1 \\ L ,, 

when 1 < p < 2. In both cases, we used Lemma [2731 Since M/N = 
o(l), we obtain that 

/ (\Vu \ p - 2 Vu - \v M Vh N \ p ~ 2 VM Vh N ) ■ Vui 
Jn 

< o{M l - n N p - 1 ) + o(l) / \Vux\ p dx. 

Jn 

Collecting these estimates together, we have proved that 

1 i$ ||5div(|?] M V/lAr| P ~VfV/;,Ar)|| LP '||VMi|| L P 

+ o(M 1 ~ n N p - 1 ) + o(l) [ \V Ul \ p dx. 

Jn 

We claim that as M — > oo, 

(2.4) ||5div(|7; M V/ijvrVV/iiv)||^ = o(M 1 ~ n A^ p_1 ), 

from which estimate (12.31) follows by Young's inequality. So, it remains 
to prove (12.41) . Since r\ M and h N are explicit functions, this follows 
from a direct computation. Noting that div(| V }in\ p ~ 2 V h^) = A p tiN = 
N p A p h = 0, we have 

dw(\ VM Vh N \ p - 2 VM Vh N ) = V^ 1 ) • \Vh N \ p ~ 2 Vh N 

= (p - l) V p - 2 MVv(M ■ )N p - l (\Vh\ p - 2 Vh)(N ■ ). 
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Consequently, since S(x) = x n , 

P-2, 



\Sdiv(\rjM^h N \ p - z r) M Vh N )\\ P Lpl 

<M^N P [ xt I \Vh{Nx)\ p dx 
JB{o,i/M)nn 



oo 



<M~ N 1 " 1 -— / / xr^VKNx'^Xr^dx 

JO J\x'\<l/M 

KM^N 1 " 1 -^ [ [ xf 1 : e~ pXn dx' dx n 

JO J\x'\<l/M 

= 0(M^ T ~ n+1 N p ~ 1 ~^ r[ ). 

This is 0{M 1 - n N p - 1 {M/N)^) = o(M 1 ~ n A^ p - 1 ) as required. This 
finishes the proof of (12 .4p . and hence that of (12. 3p . Now the lemma 
follows easily from (12. 3p . When p > 2, we have 

1= I (|Vu| + \Vuq\) p ~ 2 \Vu x \ 2 dx > [ iVu^dx, 
Jq Jq 

which, together with (12. 3p . implies the desired estimate in the lemma. 
When 1 < p < 2, we have by Holder's inequality 



3-1 
2 



{Vu^dx < n yj (\Vu\ + \Vu \) p dx 

which implies the lemma. □ 

We now prove the following result, which immediately implies Theo- 
rem [L2] in the case where the boundary is flat near the point of interest. 

Proposition 2.5. If Q is as above, there exists a sequence of explicit 
functions (%) C C^°(R n ) such that their boundary values fu = vm\bq 
satisfy supp(/ M ) C B(0, l/M) n dfl and 



lim / A 1 (f M )f M dS = 1 (0). 



Proof. If / = u \qm where u is as in (I2.ip . then (I2.2p holds true. By 
Lemma [2.31 we have 



M n-l N l-p I 7 | VMo |P rfx Cp7 (Q) 

Jq 
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where c p = p v ^ J Rn _i r](x', 0) p dx'. By ( 1A.7j) 



/ 7(|Vm| p ~ 2 Vm - |VM r _2 VM ) • Vu dx 
Jn 

< / (|Vw| + |Vm |) p_2 |Vm - Vm ||Vm | dx 
Jn 

If p > 2 then the Holder inequality and Lemmas 12.31 and 12.41 imply that 
the last expression is bounded by 

< (||Vw|| L p + ||Vw ||lp) p ~ 2 ||Vwi|| lp ||Vmo||i,p 
= o{M 1 ~ n N p - 1 ). 

If 1 < p < 2 we obtain the same estimate from 

/ (|Vu| + |Vm |) p " 2 |Vu - Vu ||Vm | dx 
Jn 

< I \Vu- Vm | p_1 |Vm | dx 



n 



^ l|VMl||rj, 1 ||Vuo 



LP II v "0|| LP 
1-n t\tp-1^ 



= o(M L - n N p - L ) 
Thus, if we define 



vm = uo 

V c P ) 

then the result follows. □ 

3. Real valued case 

Theorem 11.11 concerning the nonlinear DN map acting on real val- 
ued functions is proved in this section. When proving Theorem 11.21 we 
made use of p-harmonic complex exponentials. If p = 2 there is little 
difference between real and complex solutions, since the real and imag- 
inary parts of a complex solution can be used as real valued solutions. 
However, in the nonlinear case when p ^ 2 it is not clear how to obtain 
real valued solutions from complex ones. Therefore we cannot directly 
use the complex exponentials to establish Theorem 11.11 and the proof 
will require a real valued replacement for these functions. 

Inspecting the proof of Theorem 11.21 carefully, we see that the exact 
form of the exponential h(x) is not essential. Rather, one needs that h 
has certain properties, and it turns out that there exist real valued p- 
harmonic functions enjoying these properties. The following oscillatory 
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solutions which decay exponentially in the x n variable were used by 
Wolff [121 Section 3]. 

Lemma 3.1. The function h(x) = e~ Xn a(xi) satisfies A p h = in 
if a satisfies 

a"(x\) + V(a, a')a = 0, 

where 

v( (2p-3)(aT + (p-l)a 2 

Any solution a is smooth and periodic with period A = \(p), and satis- 
fies J a(xi) dx\ = 0. 

Note that 

\h(x)\ = e- Xn \a( Xl )l 

\Vh(x)\ = e- Xn (a( Xl ) 2 + a'( Xl ) 2 ) 1/2 . 

We denote h N (x) = h(Nx) for N > 0. 

We begin the proof of Theorem 11.11 Assume that Q C R n is a 
bounded open set with C l boundary, and xq G dQ. After a translation 
and rotation, we can assume that Xq = and the outer unit normal 
to dfl at is the nth coordinate vector — e n . We let p be a boundary 
defining function for Q, that is, p : R n — > R is a C 1 function such that 

fl = {i£l"; p(x) > 0}, 
dn = {xeR n ; p(x) = 0}. 

After scaling if necessary, we may assume that p(0) = and Vp(0) = 
e n . We define a map /:£]—»■ R™ as 

f(x) = (x, p(x)), x = (x, x n ) E Q. 

This map is invertible and close to the identity map in -6(0, 1/M) fl Q 
for M large since p is a C l function, which implies that as M — > oo 

(3.1) sup \Vp(x)— e n \ = o(l), sup |D/(a;)— J| = o(l), 
xeB(o,i/M)nn x€B(o,i/M)nn 

where I is the identity matrix. 

Similarly as in Section [2], we employ an approximate p-harmonic 
function 

u (x) = r] M (x)h N (x) 
where f]M{x) = rj(Mx) is the same cutoff function as in Section 2, and 

h N (x) = h N (f(x)) = e- Np ^a{N Xl ), 
where a is as in Lemma [3.11 and is not identically zero. 
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Note that uq G C^.{R n ). We choose a solution u by requiring that 

div(7(x)|Vw| p ~ 2 Vw) =0 in ft, 
u = Uq on dVt. 

Writing / = Mo|an, we observe that the following real valued analogue 
of Q22J holds: 

(3.2) / A 7 (f)fdS= [ 7 |V Mo | p ^ 
Jan Jn 

+ [ 7(|Vw| p - 2 Vw - |Vm | p ~ 2 Vm ) • Vu dx. 
Jn 

As in Section 2, will recover the value of 7 at by taking the limit of 
this identity as M — > 00. To estimate the limits, we need the following 
lemmas, Lemma 13.21 and Lemma 13.31 which are analogous to Lemma 
12.21 and Lemma 12.31 respectively. As in Section 2, M/N = o(l) as 
M — > 00 and 5(x) = dist(x, d£l). 

Lemma 3.2. Let C G C~(5(0,2)) and /3 > 0. Let # : R -»■ R fee a 

bounded periodic function with period A > 0. TTien as M — > 00 

M n_1 iV / C(Mx)e" p ^ (:!; ^(Arx 1 )rfx^ - I ((x',0)dx', 
Jn P ii™- 1 

5(xf((Mx)e- pNp(x) dx = 0(M 1 ~ n N~ 1 ~^), 



where c = j J g(i) dt. 

Proof. We make the change of variables x = f~ l (y), and (13. ip together 
with easy computations shows that it is enough to prove these claims 
when Q is flat near 0. We only prove the first claim and compute 

M n-i N I ^My)e- pNvn g{N yi )dy 



C(y\^yn)9(^yi)dy') e~^dy n 
\J\y'\<2 iv M J 

00 / r M N \ 

(C(y', 0) + 0(-y n ))g(- yi ) dy' e~™ dy r , 

y'\<2 iV M J 
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Writing R = N/M and ip(yi) = J Rn _ 2 ((y h y", 0) dy", we have 

1 f 00 

C(y', O)g(Ryi) dy' = - ip(t/R)g(t) dt 

E™- 1 -ft J-oo 

i 00 /.(fe+l)A 

- / ^{t/R)g{t)dt 

a' 



A 



£ ^((t + kX)/R) B 

,k=—oo 



g(t) dt 



where the expression in brackets is a Riemann sum converging uni- 
formly to L„_i C(y') 0) rfy' as M — > oo. The result follows. □ 

Lemma 3.3. We have as M — > oo 



M n-i N i- P I 7 |Vu | p rfa;^c P 7(0) 



where 



c p = — [ r]{x', Of dx', K = \ t (a(t) 2 + a\tf) 1 dt. 

V J-Rn-l A Jq 

We also have 

[ \Vu \ p dx = O(M 1 - n N p - 1 ), 
Jn 

[ \VMVh N \* 'dx = 0{M l ~ n N p - 1 ), 
Jn 

[ |Vw - r] M Vh N \ p dx = 0{M l - n N p -\M/Nf). 
Jn 

Proof. Noting that V«o = MVi/(M ')tiN + ^mV/itv, the last three 
estimates follow from Lemma 13.21 For the first claim we observe that 

Vh N (x) - Vh N (f(x)) = -Ne- Np{x) (y P (x) - Vp(0))a(N Xl ). 

Thus by the inequality (IA.60 . when computing the limit of f n j\ Vu \ p dx 
as M — > oo we may replace 7 by 7(0) and V«o by t)m( • )V/i/v(/( ■ )). 
Lemma [3.21 gives the required expression for the limit. □ 

Next we write u = u + U\ and show that Viti is asymptotically 
smaller than Vwo in the L p norm. 

Lemma 3.4. As M — > 00 



/ iVu^dx = o{M l - n N p - 1 ). 
Jn 
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Proof. The proof is analogous to the proof of Lemma 12.41 We start 
with 



/= [ (\Vu\ + \Vu \) p - 2 \V Ul \ 2 dx 
Jn 



< / 7 (|V M | + |Vu |) p - 2 |V Ml | 2 cfe 
Jn 

< [ 7 (|Vw| p ~ 2 Vm- \Vu \ p - 2 Vu ) ■ (Vm- Vu )dx 
Jn 

= - 7|Vw r" 2 Vw • Vu x dx 
Jn 

= -7(0) / I Vu | p - 2 Vu ■ V«! dx + o((M 1 -W 
in 



LP- 



In the last step we used Lemma 13.31 the continuity of 7 and the fact 
that m is supported in B(0, 1/M). 
We continue by writing 



Viio| p 2 Vuq -Vuidx = / \r] M 'Vh N \ p 2 r] M 'Vh N ■ V«i dx 



+ / (\Vu \ p - 2 Vu - \vM^h N \ p - 2 7] M Vh N ) ■ Vuidx. 
Jn 



For the last term, we can apply Lemma 13.31 and proceed in the same 
way as in the proof of Lemma 12.41 to obtain that 



v« r~ 2 v 



uo - \f]M^h N \ p 2 rj M Vh N ) ■ Vmt dx 



= o((M 1 -"^- 1 ) J T i )||Vu 1 || LP . 

For the first term in the right hand, we estimate in the following way. 
Note that 1in(x) = /i at (/(#)). We make a change of variables and use 
that V/itv(x) can be replaced by V/iat(/(x)) up to small error to obtain 



\i]M^h N \ p 2 r] M 'Vh N ■ Vmi dx 

v P M\r\ymh N \ p - 2 Vh N -Vu 1 U-\y))dy 

+ o((M 1 -"A^- 1 ) £ ? i )||Vu 1 || L P. 



Here we also used (13. ip . Then we can estimate the first term in the 
right hand side by Hardy's inequality in the same way as in the proof 
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of Lemma I2.4L Eventually we obtain the required estimate for /, 

I = o((M 1 ~ n iV p - 1 )^ 1 )||Vwa|| L p = o{M l - n N p ~ l ) + o(l) / |V«i| p dx. 

Jq 

We can conclude the proof as in Lemma [2.41 □ 

Theorem 1 1 . 1 1 now follows from the next result. 

Proposition 3.5. If Q is as above, there exists a sequence of explicit 
real valued functions (%) C C^iW 1 ) such that their boundary values 
fu = v M \ d u satisfy supp(/ M ) Q B(0, 1/M) n dtt and 



lim / A 7 (f M )f M dS = j(0). 
4 ^°° Jdn 



M^oo 

Proof. As in the proof of Proposition 12. 5[ we set 

vm = u , 

V c p J 

where the constant c p is the same as in Lemma 13.31 Then the result 
follows. □ 



Appendix A. Inequalities and wellposedness 

We first record some standard inequalities for pth powers of vectors. 
The first set of inequalities is valid for 1 < p < oo and for all £, rj G M n . 

(A.i) \v\ p >\t\ p +p\tr 2 £-(v-Z), 

(A.2) - \r,\'\ < P (\cr i + ivr'm -vi 

(a.3) \\tr 2 t-\vr 2 v\<m + \v\r 2 \t-ri\, 
(a.4) (iei p - 2 e - \vr 2 v) ■ (e - v) ~ del + Mri£ - v\ 2 - 

These results easily imply similar inequalities for complex vectors upon 
separating real and imaginary parts. If z, w £ C™ and 1 < p < oo we 
have 

(A.5) \w\ p > \z\ p + p\z\ p ~ 2 Re[z ■ (w - z)], 

(A.6) \\z\ p - \w\ p \ < pdz^- 1 + \w\ p ~ l )\z - w\, 

(A.7) \\z\ p ~ 2 z - \w\ p - 2 w\ < (\z\ + \w\) p - 2 \z - w\. 

(A.8) (\z\ + \w\) p - 2 \z - w\ 2 ~ Re [{\z\ p ~ 2 z - \w\ p - 2 w) ■ {z - w 
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Assume that ft C M n is a bounded open set, and that 7 G L°°(ft) 
is a positive function. Let IK be either R or C. We consider the well- 
posedness of the Dirichlet problem for the p-Laplace equation, 

div(7(x)|Vw| p - 2 Vw) =0 in ft, 
u = f on <9ft. 

We look for a weak solution u G W^ 1,p (ft; K), so that 

/ 7 1 Vu\ p - 2 Vu -Vipdx = for all cp e W^ft; K). 
Jn 

Here / G W 1,p (ft; K), and the boundary condition is interpreted so that 
u — fE W ' p (ft; K). We think of / belonging to the abstract trace space 
W lj? (n\ K)/Wo' p (n; K), and also write u| sn = /■ If ft has sufficiently 
nice (say Lipschitz) boundary, the trace space can be identified with 
the Besov space [23] 

W 1 ' p (ft;K)/W 1 ' p (ft;K) « S^ 1 ^ (9ft; K). 

Proposition A.l. Given any f G W l,p (Q; K) t/ie above Dirichlet prob- 
lem has a unique solution u G H /1,p (ft; K) satisfying 

wzt/i C independent of f . 

Proof. The proof is completely standard and is usually given for the 
case IK = M. [30J . We show that the same proof works for IK = C. 

We will show the solution is obtained as the unique minimizer of the 
energy functional 

E{v) = [ -f\Vv\ p dx 
Jn 

among all v in the set A = {v G W 1 *^; C);v-fe W 1,p (ft; C)}. In 
fact, let u be a minimizer. Then for any ip G Wq ' p (ft; C) 

E(u + tcp) = [ ~?(\Vu + tV(p\ 2 ) p/2 dx 
Jn 

= [ 7 (|Vw| 2 + 2tRe(Vw Vp) + t 2 \V<p\ 2 ) p/2 dx 
Jn 



and 



d f 

—E(u + t<p) =p 7|Vw| p-2 Re(Vw • V^) dx. 
dt t=o J n 



Since u is a minimizer the last expression is zero, and choosing ip purely 
real or purely imaginary gives that u is a solution, 

div( 7 |Vu| p - 2 Vu) = 0. 
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Conversely, if u is a solution with u — / G Wq ,p (£1; C) then for any 
v & Awe have by (IA.5j) 

= / 7 |Vt;| p da; 
Jq 

> [ (\Vu\ p + p\Vu\ p - 2 Re[Vu- (Vv -Vu)}) dx. 
Jq 

Since v — u G V^ 1,p (f2; C) one has i£(i> ) > -E'(m) and u is a minimizer. 

It is thus enough to show that the energy functional has a unique 
minimizer in A. We begin with uniqueness, and suppose that Ui,v,2 G 
*4 are two distinct minimizers. Then Vtti 7^ Vu 2 in a set of positive 
measure (since otherwise U\ — U2 G Wq ' p (f2; C) is constant on each 
component and thus zero by the boundary condition). The function 
Ml +" 2 also belongs to A, and the strict convexity of x h-> |x| p implies 
that 



< 



|Vmi| p + |Vu 2 | 



where the inequality is strict in a set of positive measure. We obtain 



< 



t|v«i| : 



dx < 



1 



VU\ + V"U 2 



dx 



7- 



|Vm 2 | p 



dx. 



This implies E(ui) < E(v,2) < E(u\), which is a contradiction. 

To show the existence of a minimizer, we let (vj) C ^4 be a sequence 
such that -E'(fj) — » -Eo where _E = inf ue _4 E'(u). Since / G ^4, we have 
Eq < 00 and 

/ I Vfj| < C < 00. 
in 

Using that Vj — f G W ' p (f2; C), the Poincare inequality implies 

||% - /IUl, < C\\V( Vj - f)\\ LP <C<oo. 

By weak compactness there exists u G W 1,p (ft; C) with t> 3 - — ^ u and 
Vfj — ^ Vm weakly in L p (f2; C), and since W ,p is closed we have u—f G 
W 1,P (Q;C). The function w is a minimizer: by (1A.5|) 



£K0 



7 1 Vf ,-| p dx 



> 



[ 7 (I Vu\ p + p| Vm| p " 2 Re [Vu • (Vv 3 - - V«)]) dx 
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and by weak convergence E = linx^oo E(vj) > E{u). The norm bound 
for u follows from the fact that E{u) < E(f) and from the Poincare 
inequality. □ 

Finally, we discuss a nonlinear Dirichlet-to-Neumann map (DN map) 
for the equation considered above. Write X = W 1,p (Q)/Wq' p (Q) and 
X' for the dual of X. If / e X, the DN map is formally defined by 

A 7 (/) =i\Vu f \ p - 2 Vu f -iy\ dn 

where Uf e W l,p (Vi) is the unique solution of div(7| Vu\ p ~ 2 Vu) = in 
fl with u-fe Wq' p (Q), and v is the outer unit normal vector of dVt. 
A formal integration by parts gives that 

(A 7 (/),<?>an= / A^(f)gdS= [ ^\Vu f \ p - 2 Vu f ■ Vg dx. 
Jan Jn 

The last identity can be taken as the weak definition of the DN map 
for / e X. One has 

\(A,(f),g) 9 n\ < CIIV^II^HV^IU, < C\\f\\ p w l\\g\\w^. 
Consequently A 7 maps X to X' and 

\\^U)\\x><c\\f\\ p -\ 

Some earlier results for Calderon type inverse problems for nonlinear 
equations have been based on studying the Gateaux derivatives of the 
nonlinear map A 7 at boundary values f — z where z is a constant. 
However, for the equation above the solution with boundary values 
z + tf (where t > 0) is u z+t f — z + tuf, and 

A 7 (^ + */) = 7 |V(z + tu f )\ p - 2 V(z + tu f ) ■ u\ dn = t p - 1 A 7 (/). 

Consequently the Gateaux derivatives of A 7 at z do not even exist if 
1 < p < 2, and also when p > 2 all the higher order Gateaux derivatives 
which exist are either or equal to A 7 . Thus for the p-Laplace type 
equation considered in this paper, arguments based on linearizing the 
map A 7 at constants do not yield any information and one needs other 
methods. 
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